Asymptotic solutions to the Gross-Pitaevskii gain equation: 
growth of a Bose-Einstein condensate 



O 

o 

(N 

o 
O 



I 

o 
o 



in 
o 



o 



X3 
O 

o 



X 



p. D. Drummond and K. V. Kheruntsyan 
Department of Physics, The University of Queensland, 
Brisbane, Queensland 4072, Australia 
(July 27, 2000) 



We give an asymptotic analytic solution for the generic 
atom-laser system with gain in a D-dimensional trap, and 
show that this has a non- Thomas-Fermi behavior. The effect 
is due to Bose-enhanced condensate growth, which creates 
a local density maximum and a corresponding outward mo- 
mentum component. In addition, the solution predicts am- 
plified center-of-mass oscillations, leading to enhanced center- 
of-mass temperature. 

PACS numbers: 03.75.Fi, 05.30.Jp, 03.65.Ge. 



I. INTRODUCTION 

The description of Bose-Einstein condensate (BEC) 
growth has become important due to the need to 

understand the physics of atom lasers . These recently 
developed devices which emit coherent wave-like beams 
of atoms promise a new generation of precision measure- 
ments, applications in nanotechnology and novel tests of 
fundamental concepts in quantum theory. While it is 
possible to perform first-principles numerical simulations 
of the relevant equations ||^ , a great deal of physical in- 
sight can be obtained from an analytic solution. 

In this paper, we give an analytic asymptotic so- 
lution to the Gross-Pitaevskii equation describing the 
early stages of condensate growth in a trap. The phys- 
ical insight we obtain from this is that a growing non- 
equilibrium condensate has a non-uniform momentum 
distribution across the condensed region. As a result, the 
observed density behaves as though the trap frequency is 
increased, relative to the usual Thomas-Fermi behavior 
of an equilibrium BEC. In addition, our solutions show 
center-of-mass oscillations whose kinetic energy is ampli- 
fied with the condensate growth. This implies that, while 
BEC's formed through evaporative cooling may have a 
low temperature for their internal degrees of freedom, 
the temperature for the center-of-mass motion is higher, 
leading to increased noise in the velocity and direction of 
atom laser beams. 



II. GROSS-PITAEVSKII GAIN MODEL 

We start by considering a commonly used model of a 
one-component trapped Bose-Einstein condensate - the 
Gross-Pitaevskii (GP) equation Q modified by a linear 
gain term g [b|, of the form: 



9^'(x,t) 
dt 



' Zm 



(1) 



Here ^(x, t) is the mean- field amplitude (so that 
|4'(x, i)p is the particle number density), m is the atomic 
mass, and U is the effective interaction potential. In the 
treatment of D = 1, 2 or 3 space dimensions, U is given 
hy U — AnhaL^^^ /m, where a is the scattering length 
and L is the confinement length. The potential term 
y(x) is due to an optical or magnetic trap, which we as- 
sume is harmonic. In the simplest rotationally symmetric 
case, the trap potential is given by ^(x) — mcj^|xp/(2fi,), 
where to is the trap oscillation frequency. 

It is helpful to understand the physical implications 
of Eq. (|^) by an analysis of the relevant time scales. 
In the present case, there are three different relevant 
time-scales. These are: the time-scale for the condensate 
growth, tg — 1/g, the trap oscillation time-scale, = 
l/uj, and the "healing" time-scale, th — l/(C/|\I'(xo, t)p), 
associated with the mean-field interaction potential at 
the condensate center-of-mass Xq. 

Since current BEC's typically have a relatively high 
density, the healing time is usually the smallest, which 
results in a Thomas-Fermi (TF) type equilibrium con- 
densate. If the gain time-scale is long enough, then one 
may expect that even during growth, the condensate will 
adiabatically follow the TF solution. This is commonly 
assumed in analyzing experimental data j^]. Another 
possibility is that the gain time-scale is shorter than the 
trap oscillation period, in which case we should no longer 
expect adiabatic TF like behavior. 



A. Asymptotically growing solution 

To show this distinction, we now proceed to give an an- 
alytic non-equilibrium solution of the GP equation with 
gain. Firstly, we expand the field in terms of the ampli- 
tude and phase. 



*(x,t) = A(x,t)e-^'^("^*VVC7, 
and obtain the following coupled equations: 
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dt 
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Next, we wish to investigate possible asymptotic so- 
lutions for long times - i.e., steadily growing solutions, 
valid some time after initial nucleation of the conden- 
sate, yet before any gain saturation has occurred. We 
consider first the rotationally symmetric case, where 
= + [{D - l)/r]dr, with r = |x|. Following the 
construction successfully used recently in optical fiber 
environments [p^ - and extending these to the current 
multi-dimensional case of a trapped BEC - we suppose 
that the amplitude has a self-similar behavior at large 
time. The phase is assumed to depend on the atomic 
density at the condensate center, and to have a uniform 
chirp giving a radially dependent outward momentum. 
Thus, we seek a solution in the form of: 



A^\{t)f{r/\{t)). 



2~gT 



■mg 2 



(5) 
(6) 



Here A(t) is a scaling function, while r and g are unknown 
coefRcients. Using this construction, we find, from Eqs. 
(a and (0): 



dXit) 



dt V A 



.gA/ 



(7) 



where we have defined uj'^ = uj'^ + g'^ . 

From the first (amplitude) equation, the solution re- 
quires the twin conditions that: 



(9) 



This implies that the scaling function A grows exponen- 
tially with time, having a solution of A(i) — exp(^i), 
where we can immediately solve for the growth rate, since 
clearly g = 2g/{D + 2). The physical interpretation of 
this equation is rather straightforward; the growth of the 
amplitude at any radial point is reduced below the single- 
mode growth rate g, due to an outward flow of atoms, 
which transfers part of the increased density to a larger 
radius. This effect increases with dimensionality of the 
space. 

The phase equation can now be simplified using the 
fact that the last two terms involving derivatives have 
terms in 1/A, and hence become exponentially smaller 
than the earlier terms, at long times, as the healing time 
becomes smaller. This immediately gives the following 
solution for /: 



/ = 



v/(l - |x|Vi?(i)2)/r, for |x| < R{t), 
0, for |x| > R{t), 



(10) 



where R{t) is the maximum radius given by: 
R{t) = eSV2fi/(wr(I;2). 



(11) 



The remaining unknown constant r is obtained by eval- 
uating the integral N(t) = J |*(x, i)pd^x for the total 
number of particles: 



N{t) 



2ft 



D/2 



2Cdt 



-{D+2)/2 



UD{D + 2) 



(12) 



Here Cd is a dimension-dependent constant, with Ci — 
2, C2 — 27r, and C3 = 4tt. This gives the expected 
result that the overall atom number grows exponen- 
tially with a gain constant of 2g — g{D -I- 2); that is, 
N{t) — 7V(0) exp(2gt). The reason for this apparent 
difference in gain constants is that the relatively slower 
growth in density at the center of the condensate is ex- 
actly compensated for by the increase in condensate ra- 
dius with time. 

The constant t is solved in terms of iV(0), so that: 



D/2 



2Ct 



UD{D + 2)N{0)CjD 
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(13) 



This coincides with the healing time for a condensate 
with atom number A'^(O), in TF equilibrium - but with 
an effective trap oscillation frequency Cj instead of lo. 
The final result for the particle number density 

\^{^,t)\^ ^xyyu is: 



l*(-'^)l =7^(1- 



R{ty 



(14) 



While this gives an asymptotic solution in terms of the 
initial atom number, the solution needs to be compared 
with the usual TF solution to understand the physical 
implications. In particular, we can notice an interesting 
scaling behavior for the radius, as a function of the trap 
frequency. In the present case, the radius scales as: 



i? CX (:>-2/(D+2) ^ ^^2^~2yl/(D+2)^ 



(15) 



In equilibrium TF solutions, there is a similar behav- 
ior, except that there is no outward momentum term. As 
a consequence, the TF radius {Rtf oc aj~^/(^+^)) is al- 
ways larger at a given total atom number, and the density 
is lower at the center, than in these non-equilibrium so- 
lutions with gain present. The physical reason for this is 
simply the rapid Bose-enhanced increase in number den- 
sity (and hence pressure) at the condensate center during 
the nucleation process. This effect becomes appreciable 
when g >> uj, 01 tg << 1/uj, so that the time-scale for 
condensate growth is faster than the trap oscillation pe- 
riod. 



B. Numerical results 

We now turn to comparisons between the analytic 
asymptotic solution and exact numerical results. In Fig. 



2 



1 we plot the r.m.s. radius f{t) of a growing BEC versus 
t. The result found from the above asymptotic solution 
in 3-D is simply given by f{t) — y^3/7R{t). This is rep- 
resented by the dashed line, for g = lOu! and N{0) = 10. 
For comparison, the TF r.m.s. radius for the same values 
of N{t) is represented by the dotted line. The full lines 
correspond to the results of direct numerical simulations 
of the GP equation with gain, with initial Gaussian wave 
functions of different widths. The values of other param- 
eters are chosen to correspond to a *^Rb BEC (m = 1.44 
kg, t/3 = 5 X 10"^'^ m^/s) in a trap with u)/2'k = 100 Hz. 
As one can see from the graph, despite the initial dif- 
ferences, the mean radii approach the same asymptotic 
value of ^ y^3/7R{t), which is different to the TF result 
for the same total number of particles. 

The outward momentum density p = |p| (where p = 
-i;i[^'*(V*)-iJ.c.]/2) is plotted in Fig. 2, showing that 
in a growing BEC the flow of atoms from the trap centre 
increases initially with the distance, reaches a maximum, 
and vanishes at r > R. 



III. ASYMMETRIC CASE AND 
CENTER-OF-MASS OSCILLATIONS 

The above rotationally symmetric results can easily be 
generalized to asymmetric cases, where the trap oscilla- 
tion frequencies LUi in each space direction are different. 
In this case, the trap potential term in Eqs. and (Q) 
is replaced by 



2 2 



(16) 



Assuming again that the amplitude has a self-similar 
behavior at long times, one can use the previous con- 
struction. In addition, we find that a further more gen- 
eral solution can be obtained. Following the approach 
of Kohn's theorem [|ll[, generalized to include gain, we 
allow for center-of-mass oscillation of the growing con- 
densate, independent of the inter-particle interactions. 

Thus, we seek a solution in the form of 



A = mf (Ax(t)/A(t)) 



2gT 



Ax(i) •Xo(i) 



where Ax{t) = K — xo{t), xo(i) is the center-of-mass, and 
the dot stands for the time derivative. 

Using the same procedure as before, we first find - 
from the amplitude equation - that the function X{t) is 
given by 



X(t) = exp(5i). 



(17) 



where 5 = 2g/(i) + 2), as previously. 

To treat the phase equation, in which we neglect the 
last term (x 1/A that becomes exponentially small at long 
times, we assume that 

f=Fi+F2{Ax)\ 



where the functions Fi and F2 are to be found by equat- 
ing the terms in powers of Ax^. From the terms in Axi, 
we find that each component of xo(t) satisfies equation 
XQ,i + tOiXo^i — 0, so that the condensate center-of-mass 
oscillates according to 

xoA^) = ao,i cos(wit -t- (5o,j), (18) 

where ao,i and 5o,i are the initial amplitude and phase. 

By equating the terms in (Axi)'^ and the zeroth-order 
terms, respectively, we obtain an equation for 0o(Oi 

• m v-^"^ r/ \2 2 1 

as well as solutions to Fi and F2 ' 
Fi = l/r. 



Fo = 



2hX(t) 



D {ojj + ~g^)Ax,{t) 



{Ax{t)f 



Combining these together and using the solutions for 
the center-of-mass coordinates, we finally obtain that the 
solution for the function 0o(i) is given by 



00 (t) = -(m/2ft)(xo(i) -9x0(0), 
while the solution for / is 



R^{tf 



(19) 



(20) 



in the region of space where the expression in large brack- 
ets is positive, and / = outside that region. 

The maximum radius Ri{t) in each space direction is 



R^{t)=ea' 



where we have introduced 1 



2h/imTuf) 



(21) 



I ujf +g'^. In addition, the 
constant r is solved, as previously, in terms of A^(0) by 
evaluating the integral for the total number of particles 
N{t) = N (0) exp{2gt) . The resulting expression for r is 
given by the same equation as before, Eq. (13), except 



that Lj^ is replaced by IljLi ^j- This again corresponds 
to the healing time in a TF condensate, with A^(0) par- 
ticles and effective trap oscillation frequencies uji. 

The final result for the particle number density is now: 



l*(x,t)P-^fl-E^^^^3^). (22) 



tU 



An example showing growth of the BEC, while the con- 
densate center-of-mass oscillates, is shown in Fig. 3. 
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IV. SUMMARY 

The physical interpretation of these results is that the 
asymptotic solution with gain has a density distribution 
that is similar to the TF solution, except that the trap 
oscillation effective frequency is increased, with uif — > 
ojf = (jjf + g^. For a given gain constant, this has the 
strongest effect for weakly trapped (low-frequency) direc- 
tions, in an asymmetric trap. In addition, the solution 
gives an outward momentum component due to bosonic 
stimulation effect which is stronger in the center of the 
trap. 

In addition, the results obtained show that Bose stim- 
ulation can occur to moving condensates as easily as to 
stationary ones. In cases like this, the amplitude of the 
center-of-mass oscillation does not change during the con- 
densate growth, while the total mass of the condensate 
increases. This means that the condensate center-of-mass 
kinetic energy, Ek — N (1)171 J^i^f'^o ii can grow expo- 
nentially to large values - even in the absence of tech- 
nical noise. While the present model of gain is simpli- 
fied, similar types of center-of-mass motion are found in 
first-principles simulations of condensate formation 
This suggests that, while EEC's formed through stimu- 
lated emission may have a low temperature relative to the 
center-of-mass, the center-of-mass motion may itself have 
a higher temperature. Current experimental measure- 
ments of single atom velocity distributions | [T^ appear 
insensitive to the center-of-mass temperature of the con- 
densate. The effect is analogous to the increased uncer- 
tainty in the frequency and pointing stability of a multi- 
mode optical laser, compared to a single-mode laser. 

We emphasize that the asymptotic solutions given here 
are only applicable for condensate nuclei that have al- 
ready formed, as the spontaneous emission noise is omit- 
ted. The gain medium is assumed to be unsaturated 
and to equilibrate rapidly, giving a uniform gain con- 
stant across the growing condensate. In addition, one can 
expect damping to occur due to interactions with non- 
condensed atoms which are not treated here. These 
interactions, however, can only equilibrate the temper- 
atures of condensed and non-condensed center-of-mass 
motion, rather than providing additional cooling of the 
condensate center-of-mass motion. Other damping mech- 
anisms that will intrinsically be present in a BEC are 
two-body losses, which effectively change g to 5 — Fj^'p, 
where F is the two-body loss rate. While these affect the 
form of the asymptotic solution, two-body (or higher- 
order) losses do not change the amplitude of the conden- 
sate center-of-mass oscillation 

In summary, we have found an asymptotic solution to 
the Gross-Pitaevskii equation with gain, which has the 
advantage of yielding an explicit analytic result of great 
physical transparency. The solution shows that the non- 
equilibrium behavior of a growing Bose-Einstein conden- 
sate generally includes an outward momentum compo- 
nent and spatial oscillations. As a result, we suggest that 
the state of a trapped BEC is not a canonical ensemble, 
and should be characterized by at least two distinct tem- 



peratures: one for the internal degrees of freedom, and 
one for the center-of-mass motion. 

The authors acknowledge the Australian Research 
Council for the support of this work. 
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FIG. 1. The r.m.s. radius f (in a logarithmic scale) of 
a growing condensate versus t, for g = Wuj. The full lines 
(1), (2) and (3) correspond, respectively, to initial Gaus- 
sian wave functions having r.m.s. radii tliat are twice larger, 
equal, and twice smaller than the asymptotic r.m.s. radius 
f(0) = Y^3/7ii(0). The dashed line is the asymptotic result, 
while the dotted line corresponds to the equilibrium (TF) so- 
lution. 




r(nm) 

FIG. 2. The momentum density p versus the radial dis- 
tance from the centre of the trap r. The full line corresponds 
to the resulting distribution of case (1) in Fig. 1; the dashed 
line is the corresponding asymptotic solution. 
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FIG. 3. Growth of a BEG with a center-of-niass oscillation 
present. Shown is the condensate density |^'(a;,i)p, as found 
from the asymptotic solution in a symmetric trap, with a cen- 
ter-of-niass oscillation in the x-dircction and N(0) — 10. The 
gain coefficient is chosen as g = O.lu), while ao,i = 4.86 yum, 
and So,i = iv/2. Other parameter values are as previously. 
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